arXiv:1508.04844vl [math-ph] 20 Aug 2015 


Euler Polynomials and Identities for Non-Commutative Operators 

Valerio De Angelis and Christophe Vignat 
Department of Mathematics, Xavier University of Louisiana, New Orleans, 
Department of Mathematics, Tulane University, New Orleans, 

L.S.S. Supelee Universite Paris Sud Orsay 

Three kinds of identities involving non-commutating operators and Euler and Bernoulli polyno¬ 
mials are studied. The first identity, as given in 0 by Bender and Bettencourt, expresses the nested 
commutator of the Hamiltonian and momentum operators as the commutator of the momentum 
and the shifted Euler polynomial of the Hamiltonian. The second one, due to J.-C. Pain Q, links 
the commutators and anti-commutators of the monomials of the position and momentum operators. 
The third appears in a work by Eiguieira de Morisson and Ering [Gj in the context of non-Hermitian 
Hamiltonian systems. In each case, we provide several proofs and extensions of these identities that 
highlight the role of Euler and Bernoulli polynomials. 
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I. INTRODUCTION 


Special functions appear as a natural tool in many areas of theoretical physics; in quantum physics, 
Hermite polynomials are the natural basis to describe the behavior of the quantum harmonic oscillator. 
The Euler and Bernoulli polynomials are ubiquitous in number theory and combinatorics; in physics, 
they appear in various ways, such as in their association with zeta functions and in their rich interplay 
with particle physics; see [3 for examples such as the Casimir effect and string theory. 

Bernoulli and Euler polynomials also appear in the field of non-commutative operators in quantum 
physics, in more subtle ways; for example, the Bernoulli numbers appear in the the linear term of 
the celebrated Baker-Campbell-Hausdorff formula. Their role in quantum algebras is also detailed in 
0 - This paper describes three distinct contexts of quantum physics in which the Bernoulli and Euler 
polynomials play an important role. 

Define the commutator of two operators p and q as 

[p, q]=pq- qp 

and their anti-commutator as 


{p, q} =pq + qp. 

Define moreover the nested (or iterated) commutators 

[p, q ]2 = [[p> 9]: 9]: [p, q]3 = [[[p> 9], 9], 9] = [[p, 9)2:9] 

and more generally 


{p^q]n = [[p>9]„-i ,9] , n>l 


and accordingly the nested anti-commutators as 

{P. 9)2 = {{P. 9} , 9} , {P, 9}3 = {{{P: 9} , 9} , 9} = {{P, 9)2 > 9} 

and more generally 


{P, q}n = {{P> 9}„_i, 9} , n>l. 


For example 


{p, 9)2 = {P9 + 9P, 9} = P9^ + ‘^qpq + q^p 


and 


[P> 9)2 = [P9 - 9P, 9] = Pq^ - ‘2qpq + 9^P- 


It can be checked by induction on n that the general cases are 



and 



as a consequence, the exponential generating functions of these two sequences are 



n>0 


n>0 
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Note that [p, q] = [p, q]^ and {p, q} = {p, q}^. 

In the sequel, we will consider the more particular case where the non commutative operators p and 
q satisfy the identity 

[q,p]=qp-pq = i (I.l) 

where = — 1. A more rigorous notation would be 

qp — pq = il 


where / is the identity operator, but in the absence of ambiguity, we will systematically use the 
simplihed notation ini. For example, take q = and p = x in the sense that q and p act on any 
differentiable function / as follows 


and 


Qf {x) 



*/' (x) 


pf (x) = xf (x). 

Then {qp — pq) f = if so that indeed 


qp-pq = 

in the sense that, for any differentiable function /, 

{qp -pq)f = if. 

However, we wish to prove the results stated in the next sections in the general case of two arbitrary 
operators p and q, assuming only that p and q satisfy jni). 

Let us introduce moreover another operator, the Hamiltonian operator id, as 

H = \{p' + ■ ( 1 - 2 ) 

Our aim in this paper is to study three types of identities on ^ecific commutators and anti¬ 
commutators; the first one was studied by Bender and Bettencourt [3, and links the n—nested anti¬ 
commutator of the operators q and H to a simpler anti-commutator of q and a polynomial version of H. 
This result was not proved, but only inferred by these authors. Only recently, a proof was provided in 
[^, based on the integral representation of operators, more precisely the generalized Weyl transform. 
Our approach consists in first transforming Bender and Bettencourt’s identity into a simpler form, 
and then in providing two different proofs, one algebraic and another based on the properties of some 
operators. A third, purely analytic proof, is obtained by realizing the considered operators and looking 
at their action on a sufficiently large set of functions. 

The second identity studied in this paper was introduced recently by Pain Q: it expresses the 
commutator of monomials in operators p and q as a linear combination of their anti-commutators. 
We show that the use of generating functions not only gives a simple proof of this result, but also 
allows to derive the converse identity, expressing the anti-commutator of monomials of p and g as a 
linear combination of their commutators. In both cases, the Euler and Bernoulli polynomials play an 
important role. 

The third identity appears in the context of non-Hermitian Hamiltonian systems, as studied by 
Figuieira de Morisson and Fring in Q: the real and imaginary part of the non-Hermitian Hamiltonian 
can be shown, under assumptions that will be detailed later, to be related by a linear identity involving 
Euler numbers. 

One of the main intents of this paper is to show that for a given identity on non-commutating 
operators, a variety of different proofs - each pertaining to a different area of mathematics - can 
be exhibited. On the way of each of these proofs, some interesting results, of either combinatorial, 
algebraic or analytical nature, may appear. 
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II. AN IDENTITY BY BENDER AND BETTENCOURT 


A. Introduction 


In 0, C. Bender and L. Bettencourt inspect the first values of the operator 


1 




for n = 0,1,2, 3,4, 5 and 6, they find respectively 


1 


1 




1 


1 

2 


. 3,5 

q,H --H + — 


1 

2 






25 

16^ 




and 


1 

2 




4 16 



This suggests the following result 

^ U’ H}n - 2 + 2 ) } 


(II.l) 


where, by a real tour de force, the authors identify En (x) as the Euler polynomial of degree n, defined 
by the exponential generating function 


E 

n>0 


(^) n 

-;-Z 


2e^"^ 

+ 1 ’ 


\z\ < 27r. 


(II.2) 


Since this identity is only inferred in 0, we propose to prove it in several ways using algebraic or 
analytic methods. 


B. Preprocessing 

As a first step of the proof, we transform the identity (III. II) into the equivalent form 

Next, we consider the following lemma. 

Lemma 1. A polynomial Pn (x) satisfies 

(x) + {x + l) = x^ (II.4) 


if and only if 


Pji i,x') — E>fi (x), 


the Euler polynomial of degree n. 
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Proof. The fact that the Euler polynomials satisfy 

En (x) + En (x + 1) = 2x", n € N, 


(II.5) 


can be deduced from the generating function pi.2|) . The fact that this is an equivalence can be checked 
as follows: starting from 

En (x) + En (x + 1) = Pn (x) + Pn (x + 1) 
which shows that (x) is of degree n, and writing 

Pn (x) = + ... 

then, taking x —>■ +oo shows that is also the leading term of En (x). Then applying the same 
argument to the new polynomial 

Pn (x) - a^n'^x^ 

shows by induction on n that both polynomials and En coincide. 

Another proof of this Lemma uses the integral representation of Euler polynomials that can be found 
as formula 24.7.6 in S, 


n + OO / ^ N » 

En {z) = J ( ~ 2 ^ {t^u) du. 


(II.6) 


Start expanding 


Pn{x) = 


/c=0 


and rewrite pL4[l as 


^ Z! ^ Z (a; + 1)'^ = 


X . 


Substitute x ^ lu — ^ + x in both sides and integrate again the function seek (ttu) , giving 

n 1 

p P Z4”^-^fc (x + 1) = En (x) . 

^ k^O ^ k^O 

Using the property pl.5[l of Euler polynomials, deduce 


= Pn (x) = En (x) 




Let US introduce next the symbolic notation, with a € 


{{q,H} + a)^ = J2 




^ — k 


U,H}k 


□ 


(11.7) 


and the convention {q, 7 J}q = g. 

We will need the following identity, the proof of which is elementary. 
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Lemma 2. For a S R and n G N, 


Proof. Take the exponential generating function on both sides. On the left-hand side 

IIS 

n>0 


while on the right-hand side 


5]({<7,iI} + a)„- = e-e^V" 

n>0 ”■ 


so that both generating functions coincide. Since they are analytic functions, the corresponding se¬ 
quences coincide. □ 

Using both previous lemmas, we can now deduce the following equivalence. 

ProDosition 3. An eauivalent form of identity ill. 1]) for am) is 


1 

2 " 


q,H 



1 

2n 




n 




Proof. Assume hrst that 

holds. Then 

where we have used the property (III. 51) of Euler polynomials. 
Reciprocally, assume that 


U,H-} 


holds. As can be checked from their generating function (III. 21) . the Euler polynomials are Appell 
polynomials, i.e. they satisfy 


hence we have 


= (O)crS 

k^O ^ ^ 

{< 7 , Er, (i?)} = I q, Q E^_k (0) I 

so that we need to compute 


/c=0 
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and to show that their sum is equal to 


q,H - 


1 


By Lemma ([2]), the first sum is 

(k) ( 0 ) ^ ({?, H} - 1 )^ . 


Applying pL7[) . this is 


A;=0 


fc=0 


i: (;;) (0) ^ E (I) («.«}.=e (';) f«. e 


Z=0 


n—l 


p—0 


-i 


En-l-p (0) 


(-If 

2P+1 


=i:; 


1=0 


f 1 


2' 


En—k I 


Thus 


and similarly 






z=o 




I 2' 


7?n.-/ —- 


k=0 

SO that the sum equals 


E(:)^-(o)i^>^n^=E 


i=0 


"tA {q^H} 


l J 2' 


^-En-i 


{q,H}i 

on / y 


1=0 


I / 2' 


Lin-; ( 9 ) + Lin-i ( 9 ) ) — 9„_1 E 


1=0 


{q.H}i ( 1 

2 


I / 2' 


n—l 


where we have used (III.5I) again. 

The last sum is now identified by Lemma ([U as 


^ {{q,H} - 1 )„ = ^ {q.H-\ 


and this is the desired result. 


□ 


We conclude that proving (IIL1|) is equivalent to be proving the following result. 

Theorem 4. If p and q are two operators that satisfy EP and H is the associated Hamiltonian 
defined as in if/.ill) . then for any integer n > 0, 

^ ({q, H} - 1)„ + ^ ({9, H} + 1)„ = {q, iL"} (11.8) 

with the symbolic notation ({g,i7} + a)^ as defined in if//. 71) . 

Example 5. The case n = 1 reads 

2 ({q^^} “ l)i + 2 ({?,-ff} + l)i = {gjL/} . 

This is easily checked by expanding the left-hand side as 

t {{q, //} - l)i + t ({q, iL} + l)i = i ({g, + ^ {{q, //} + !) = {g, H} . 



















Example 6. The case n = 2 reads 


^ “ 1)2 + ^ + 1)2 — . 


The left-hand side is now 

^ ({ 9 )^} “ 1)2 + ^ -I- 1)2 = - { q , H }2 + 2 { 9 : = 2 -^^2 + 9 ) ■ 

But since [qjH]^ = [[ 9 ,-ff] ,H] and [q,H] = ip, we deduce [q,H]^ = [ip,H] = i[p,H]. It then follows 
from [p,H] = —iq that [q^H ]2 = q. 

We deduce 

\ ({9, h }2 + q ) = l {U,H }2 + [q, H]^) 

= i {qH^ + 2HqH + H^q) + ^ {qH^ - 2HqH + H^q) 

= qH^ + H^q= {q,H^} 

which is the desired result. 


C. An algebraic proof 

The first proof of Thm |4] (and thus of Bender’s identity (III. Ill ') is given in the case n even only, 
replacing n by 2n in (III. 81) : its extension to the odd case follows the same lines. It starts with the 
expansion 


1 1 " 
;j({9,i?}-l)2„ + :^({9,i?} + l)2. = E 


k=0 


2n 

2k 


{q.H} 


2k 


in which is replaced by {[g', i?] 2 n- 2 fc ^^} 2 k [^’^] 2 n- 2 fc = We obtain 


/c =0 


E ( T,] = E («) E (T)"' 


2k 


fc =0 


2n 


2k 


2k 


1=0 


2k 


T2k—l 




2n 


2k 


=E(;:)E( , )"'( E 


/=0 


2k 


/'2n—2k 


r—Q 


2n — 2k 


'f' jijr ^ j:j2n—2k—r \ 


i-iy H^qH 


n 2k 2n—2k 

= EE E 

k—O /=0 r=0 

n 2k 2n 

= EE E 

k—O l—O s—r+l—0 


2n\ (2k\ f2n — 2k 
2k 


2n\ (2k\ f2n — 2k 


■■-I 


2k \ I 


{-If H^+^qH'^^-^-^ 


We thus need to show that 


n 2k 




k=0 1=0 

Differentiating the identity 


2n\ (2k\ f2n — 2k 


•.-I 


2k \ I 


(- 1 )-' = 


0 l<s< 2 n-l 

22 "-i s = 0 , s = 2n 
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i times with respect to x and j times with respect to y, and then evaluating ai x = y = 1 gives 


E 

fc =0 


2n\ / 2k\ f2n — 2k 


2k 


J 


2n\ f2n — i 

j 


c^n—i—j — 1 


if 0 < i < n, 0 < j < n and f + j < 2n — 1, while 

^2n\ (2k\ (2n — 2k\ \ (2n 


E 

fc=0 


2k 


2\n 


(1+ (-!)") 


if * = j = u. We deduce, for 1 < s < 2n — 1, 


s n 


^-EE 


= 2 


1=0 k=0 
2n-s-l 


2ri\ (2k\ f2n — 2k 


:-l 


2k \ I 


(-ir' = E(-i)' 


(2n\ f2n — I 


1=0 


■■-I 


E(-i) 


1=0 


i-z (S 


= 0 . 


A straightforward computation gives bo = & 2 n = 2^” ^. This concludes the proof. 


D. An operator-based proof 


Another version of the proof of Thm 2] is obtained now by identifying the nested commutators and 
anti-commutators and {q,H}). as powers of simple operators. 

Define the operators Ah and Bh acting on q as 

AHq=[q,H], BHq = {q,H}. 

We first verify the following identities. 

Lemma 7. The operators Ah and Bh satisfy, for k G N, 

A%q=[q,H],, BU = {q,H}, (II.9) 

and 


{AH + BH)’^q = 2\H\ {AH-BHfq = 2^H\. (II.IO) 

Moreover, Ah and Bh commute. 

Proof. The identities (III. 91) and (III.101) are easily proved by induction on k. Moreover, BnAnq = 
qH^ — H^q = AnBnq so that Ah and Bh commute. □ 


We obtain the final identity using (|ff.9|) and (III.101) as follows: 

^ ^ /c=0 ^ ^ 

= 2 + Bnf'^ + {Ah — Bh)^^^ q 

= 2"^^-^ (qH^^ + H^'^q) = 2^^^-^ {q, . 
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E. Another Proof: an analytic approach 


The third proof of Thm|3]is purely analytic: we choose the realization 


P = x, q = i—. 

ax 


This is a faithful realization [l|: any result proved with this realization holds for any other operators 
p and q that satisfy the commutation relation m- 

Next we consider the Hermite functions defined by the Rodrigues formula 


i’n (x) = (-If 


1 


\/ 2 ' 




=e 2 


n'VTT 


x2 d" 


dx" 


n e N. 


These functions satisfy the differential equation 

(x) + x'^^n (x) = (2n + 1) tpn (x) 

or equivalently 

1 


H-lpn (x) = + - j Ipn (x) , 

so that they are eigenfunctions of the Hamiltonian. Since they form a basis of i 
to verify (jll.81) on this set of functions. But it is easily verified that 


(II.ll) 

it is sufficient 


{q, (x) = Q (x) 


with, from (III.111) . 


and with 


so that 


Vi (^) = M + 


1 


Vi (it) = * \ oV'i-i (x) - 


^—k 


V’i (x) 


1 + 1 


V'i+I (x) 


{q, i?}„ (x) = Q H'^qH--% (x) 


k=0 


I - 


1 + 


n—k 


(x) - E 


k=0 


1 + 


1 + 1 


1 + 


n—k 


i’l+i (x) 


We deduce 


= (20"' i/'i-i (it) - (2^ + 2)" Oi+i (it) = *2" 2 2 (a;) _ (/ + i)"+ 2 .0^^^ ^ 

{{q, H} + 1)„ 'ipi (x) = p (pj {q,H}^ Oi (it) 

n z X 

= E (it) -{1 + 1)'=+" Oi+1 (it)) 

h—n \^/ 


= + V’i-i(it)-* 


/ + ! 


{21 + 3)”0i+i (it) 
















and accordingly 
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- 1)„ (x) = «Y ^ (2? - 1)” ^i-i (x) - {21 + 1)” V'i+i (a;) 


so that 


[({g, H} + 1)„ + ({g, H} - 1) J (x) = *- [{21 + 1)” + {21 - 1)^ (x) 


1 + 1 


[(21+ 3)"+ (21 +in V-i+i (x). 


Now 


{g, iJ"} V'/ (x) = (x) + i?”g?/'i (x) 


= '^+ 2 ' * 


^V'Z-l (x) - (x) 


+ VU'"0 V';-i(x)-*^^^(l + n V’m(x) 


which proves the result. 


III. IDENTITIES ON COMMUTATORS AND ANTI-COMMUTATORS OF MONOMIALS 

A. Introduction 


In this second part, we change notations compared to the first part: in order to follow J.-C. Pain’s 
notations, we assume now that the operators p and g satisfy the more general commutation relation 

[p, q]=pq-qp = c, 

with c S C\ {0} . Bender’s results in the former section correspond to c = —i. 

In the recent publication Q, J.-C. Pain derived the following identity between commutators and 
anti-commutators of monomials of the operators p and g. 

Theorem 8. eg. (43)]If p and q are two operators such that 

[p, q]=pq-qp = c, 

then the commutators [p",g"*] and anti-commutators {p^jg™} of monomials are related by the convo¬ 
lution identity 


(III.l) 


-pU ^m- 

E \kEk{0) 

' pU-k ^m-k 'j 

nl ’ ml 

J 

1 

{n — k)l ’(m — k)l j 


where (0) is the Euler polynomial /? / of degree k evaluated at 0, as defined by ilLB) . 
The proof given in Q relies on the fact that if the identity 


p"- q"- 
n\ ’ ml 


min(m,n) 

= E 


k'^k J P 


rji — k ^m — k 


k\ ( (n — /c)! {m — k)\ 


(III.2) 


holds for some sequence {xfe} , then these numbers Vk should satisfy (see i) 


Vk 


1^1 


+1: 


which, with vq = 1, characterizes exactly these numbers as Vk = —E^ (0) , k>\. 
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B. A proof using generating functions 

We give here another proof based on generating functions: consider the bivariate exponential gen¬ 
erating function of the left-hand side of (IIIL2I) 


E = 

n! m\ 

m,n>0 


For the right-hand side, we have the bivariate generating function 


m.in(m,n) 


E E 


m,n>0 k—1 


c'^klv. 


^ +00 


m,n>0 k—1 


k\ 


(n — fc)! ’ (m — fc)! 


^ j ^ 


k\ 


E 


\fc=l / m,n>0 

and we deduce, denoting z = cuv, 


(up)" {vqY 


n\ m\ 


^ + cx:) 


HE 


\k^l 


[cuvY 

k\ 


Vk {e“He’'«} 


^±50 


But since [p, q] = c, using McCoy’s identity [3, Formula (14)] 


k>l 


with / (p) = exp (up) and / (q) = exp (vq ), we deduce 


\^np ^vq 


] = (1 - e-) 


^np^vq 


and, since - {e“P, e"?} = [e"^, e”®] - 2e“Pe"H we deduce 

{e“P,e’’«} = (1 + 6"^) 


(III.3) 


(III.4) 


(III.5) 


(III.6) 


In fact, these two identities can be obtained without using McCoy’s identity (1111.41) as a consequence 
of the celebrated Baker-Campbell-Hausdorff formula that tells us a more precise result, namely that 


gPg? _ ^p+q+h\P’<i\ = gP+'S+f 


and 


g^gP = e«+P- = ; 


from these two formulas, we deduce 


gP+<? _ 


and then (|III.5I) and (IIIL6F 

Hence the commutators and anti-commutators are related as 

[g^P g«9i = |e“P |g“P g«n 

L J 1 _p g-^ J -f 1 ^ 

Identifying with (IIII.3I) gives 



fc=i 


(III.7) 
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or 


+°° 

fc =0 


- 1 
+ 1 


2e^ 

IT? 


+ 00 

E 


/c=0 


Ek (1) .fc 

k\ 


so that tifc = Ek (1). Using (III.5I) . we deduce, for all fc > 1, 


Vk = Ek ( 1 ) = -Ek ( 0 ). 


C. The reciprocal identity 


One of the interesting features of this generating function method is that it allows us to invert the 
identity (|III.1I) . giving an expression of the anti-commutators in terms of the commutators. First recall 
that the Bernoulli numbers are defined by the exponential generating function 


E 

n>0 


Bn 

nl 



\z\ < 2 tt. 


We can now prove the following result. 

Theorem 9. If p and q satisfy [p, q] = c then for m, n > 0, 


/ p" g™ \ 2 

pn+l qTn+l 

min(m,n) . 

1 r, C Bk+1 

pU-k 

^m — k 

[ n! ’ to! J c 

_(n -1-1)! ’ (to -I-1)!_ 

^ k\ k -\- 1 

k>l 

(n — fc)! 

(to — fc)! 


(III.8) 


where Bn is the n—th Bernoulli number. 
Proof. First transform the identity 


into the equivalent 


ecuv 2 

^CUV ^ 


- 1 
+ 1 


^up ^ I 1 + 


{e“P, 


= (! + -- " 


z — 1 


- 1 , 

e’'«] = I 1 + - E I e"’]. 


z ' k\ 

k>0 


Since Bq = 1 and Bi = — i this gives 


{e“P,e"«}= I 1 + - I 1 


k>2 ) ) 


which simplifies to 


1 B. 


{e“^e''n = 2 - + ^ 


fc+i 


z ^' fc! fc -|-1 
k>\ 


Expanding the exponentials in the commutator and anti-commutator and identifying the corresponding 
powers of u and v with z = cuv gives the desired result. □ 
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As remarked by J.-C. Pain in Q, using the fact that 

Afc(0) = -2(2'=+i-l)|^, 

the original identity (IIII.ll) can be rewritten as 


pn qm 

n\ ’ 777,1 


miii(m,n) ^ 

= 2 y: -!)«'■« 


^n—k 


r^m — k 


k + 1 \ {n — k)V {m — k)! 


which now shows more resemblance with its reciprocal version (IIII.8I) . 


D. An analytic approach 


We look now for another proof of the identity (IIII.ll) using an analytic approach: the operators p 


and q are realized as 


P = c/' (x) , qf (x) = xf (x) 


so that they satisfy the condition 


LJ^ 7 T j ''y 

and we look at the action of the commutators [p", q^] and anti-commutators {p", g™} on the monomials 
/ (x) = X* where we assume that I > n. An easy computation gives [? ] 


and 


r,n — k 


p" g"* 
n! ml 


a—k 


X = 


ml 


m +1 
n 


X 


,Z —n+m 


(n — k)\ ’ (m — fc)! 
Hence we need to prove that 


'i — k 


X = 


I — k + l 


{m — k)l\\ n — k ) \n — k 


,Z —n+m 


m +1 
n 


min(m,n 

= - ^ £ 1 + 0 ) 


m — k + I 
n — k 


or 


V n ) \n) ^ \k } W n — k ) \n — k 

equivalently, since the Euler polynomials satisfy (1) = —Ek (0) for fc > 1, 


min(m,n) 


E 


fc =0 


'm\ fm — k + I 
n — k 


min(m,n 

= E 

A;=0 


; — k 


We will prove the following more general result. 
Theorem 10. The following identity holds 


in(m,n) 

E • 

k=0 


m\ f m — k + I 
n — k 


min(m,n) 


= i: (-+!)■ 


k fm 


k=0 


k J \n — k 


As a consequence, the Euler polynomials satisfy the identity 


min(m,n) 

E 

/c=0 


m\ f m — k -\-1 


min.(m,n) 


^ \k J \ n — k J ^ 

k=o \ / \ / 

The particular case z = 0 gives the desired identity UII.SM . 


— E ^k{z + l) 


n—k 


(III.9) 


(III.IO) 


(IILll) 
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Proof. We first proceed to prove the first identity (IIII.IOI) . Defining 


a(m,n) 


pi^)= E 




m\ ( I 


min(m,n) 


E 


fc=0 


m\ f m — k + I 
n — k 


we thus need to show that P {z + 1) = Q {z). Expanding each (x + 1)^ in the expression for P {z + 1) 
gives 

min(m,n) / \ / i \ min(m,n) k 


E (=+i> 


fm\ ( I 




:j= E E 


k—0 j—0 


k\ A (m\ ( I 


j J \k J \n — k 


minfm.n) . . m — j , ✓ 

= E 7 s "‘.^1 -d- 


Considering the coefficient of in this last expression, we see that the identity to be proved is 


m- j + l\ _ I'm-jW I 
■ ^ V k )\n-k-j 


This is known as the Chu- Vandermonde identity, and is a direct consequence of the identity 

(i+tr-^+' = (i+<r-^(i+t)', 


by considering the coefficient of in each side. 

Another proof based on hypergeometric functions is as follows: identify each of the polynomials P 
and Q as a Gau£ hypergeometric function 


P(z) = 


2.?^! 


—m, —n 
I + 1 — n'' 


z] , Q{z) = 


m + I 
n 


2 F 1 


—m, —n 
—I — m 


The desired identity P {z + 1) = Q (z) is a simple consequence of the linear transformation formula 
for the Gaufi hypergeometric function that appears in as identity 15.3.6, 


2 F 1 


■z = 


r (c) r (c — a — 6) 
r (c — a) r (c — 6) 


■Pi 


a, b 

a + b — c + 1 


1 - z 


, n ^c-a-b^{c)V{a + b-c) 

+ - T ./ - 2.^1 

r (a) r (6) 


c— a,c— b 
c — a — b + 1 


1-z 


with the choice a = —m, b = —n, c = Z + 1 — n; note that the second term in the right-hand side 
vanishes since we assumed —n -f Z -f 1 > 0 and since T (—n) = -|-c». 

The last step of the proof is to transform the general identity (IIII.IOI) into identity (IIII.III) for 
Euler polynomials; this is obtained by using the integral representation pi.611 of Euler polynomials. 

Replacing z^ by [lu — ^ + z) on each side of (IIII.I0|) and integrating over u against the seek (nu) 
function gives the desired result. □ 


E. A generalization 

Exploiting the linearity of identities IlIII.ll) and (jIII.81) allows to extend them to more general func¬ 
tions. We explicit this generalization in the case of two functions / (p) and g (p) that we assume regular 
enough to have an integral representation as Fourier transforms: 

fip)= f f (u) F-^du, g{q)= f ~g (v) F^’^dq. (III. 12) 

JR JR 
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Theorem 11. Assume that the functions f and g have the Fourier integral representation Ulll.lSfl for 
the operators p and q such that [p, q] = c. Denote moreover 


F{p)= f 

Jm. 


an antiderivative of f and accordingly G (p) for an antiderivative of g. 
Then we have 


[f ip), g (g)] = -Y^E,{0)^_ {/« ip ), (g)} 


fc>l 


and 


{f{p),g{q)} = -[F{p),G{q)]+2j2 


,k r 


k>l 


Bk+i c‘ 

k + 1 k\ L' 


/W(p),g« (g) 


Additionally to McCoy’s identity Formula (If)] 


[/(p),g(g)] = -^^/Wb)5W (g) 


k>l 


we have the obvious 


{/ ip ), g (<?)} = 2/ {p) g (g) + E (P) 


k>l 


and the two last cases are 


k>0 


/W(p),g« (g)' 


and 


/(P)s(9) = 5i:^(-c)‘{/“>(p).<;“>(9)}- 

^ fc >0 


Proof. We give the proof of the first identity only, the other proofs follow the same pattern. 
The commutator [/ (p), g (g)] is computed from the integral representations (IIII.12I) as 


[[ /(M)g(r^)[i 


e“«] dudv = ff f (u) g (?;) ^ 3 ^- - e™n dudv 

i/R 2 e ^ + 1 


where we have used (IIII.7I) replacing {u,v) by {iu,iv) and remarking that 2 ; = cuv is changed to — z. 
Now writing 


e"^ -1 1 - ^ 2e^ 

e-^ + 1 “ 1 + e^ “ ~ + 1 




k>0 


k\ 


f (u) ~g (v) e“«] dudv = f (u) g {v) (l “ E 


fiu)g{v){e^’^P,F'’‘^}dudv 


- E II («) ^'5 (v) dudv. 


k>0 


we deduce 
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The first right-hand side term is simply the anti-commutator {/ (p), g (q)} . Each integral in the second 
term is identified as 


[ u'^f (u) e^^Pdu = {-tf (p) 

Jr 


and the same for g. We deduce the second right-hand side term as 


E 

k>0 


kl 


(<?)}• 


Hence the right-hand side reads 

{/ (p ),9 (<?)} - E (p) ^ ( 9 )} = - E b) ^ 9 ^’^^ (g)} • 

A:>0 ■ k>l 

The final result is obtained by remarking that i?/c(l) = (—l)^i?fc(0),fc>0. □ 


IV. NON-HERMITIAN HA MIL TONIAN SYSTEMS 
A. Figuieira de Morisson and Fring’s results 

Another context in which the Euler numbers appear naturally is the study of non-Hermitian Hamil¬ 
tonian systems: in Eiguieira de Morisson and Fring consider a ^tem ruled by an Hamiltonian H 
that is not Hermitian but PT-symmetric, hence pseudo-Hermitian |lll |: it is assumed that H is similar 
to the Hermitian Hamiltonian h, and that the similarity transformation can be realized under the form 

h = eiHe-i (IV.l) 

for some Hermitian operator q. Writing the non-Hermitian Hamiltonian 

H = ho + ihi 

where hj = ho and h\ = /ii, the pseudo-Hermitianity condition 


is expressed as 


[ho — ihi) e ^ = ho + ihi 


or equivalently as 

ho — e'^hoe~‘^ = i{hi+ . (IV.2) 

This equation is solved in hi as a function of ho in Q, looking for coefficients such that 

- Tl. 
n—1 

In order to comply with the notations in [^, note that in this section, the nested commutator in (IIV.3|) 
is now iterated from the left: 


[ho,q]n 


[q, [q,..., [q, ho]]] . 
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It is found in 0 that all even-index K 2 „ numbers vanish and that 


1 


1 


1 


17 


31 


^1=2, «^3 = -4, 2, = «^9 = y... 


(note the sign error on Kg in i). Then h is computed using (|IV.1[1 and it is found that 

Ari 




n>0 


2"n! 


[ho,q]r 


where the coefficients A„ are related to the k„ by 



These coefficients are identified as X 2 n+i = 0 and 

Ao = l, A2 = —1, A4 = 5, A6 = —61... 

from where it is concluded (without proof) that 



(IV.4) 


B. A proof and why Euler numbers appear here 

In order to prove these results, we’ll use a new method based on symbolic computation in the spirit 
of the classical umbral calculus as introduced in [T^ , see also [I^ . We define an Euler symbol £ such 
that 


£^ = Sn (0); 


the principle of umbral calculus is to perform the computations replacing the sequence of numbers 
by the powers and to obtain the final result by the reverse substitution of each f by its counterpart 
En- 

For example, 


e 


zS 


E 

n>0 



E 

n>0 


ErM^. 

n\ 


2 

e^Tl 


and more generally 


^z(,S+x) 


E 

n>0 


{£ + xf 

nl 


■y En {x) 
^ nl 

n>0 


2e^^ 

e^TT 


so that f -I- x is the symbol for the Euler polynomial, in the sense that 

(£ + x)" = En {x). 


(IV.5) 


The fundamental property of Euler symbols that explains their appearance in this context is the 
following extension of the identity (III. 51) . 

Lemma 12. For any analytic function /, 


f {z + £)+ f {z + £ + 1) — 2 f {z) ] 


(IV.6) 
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Proof. By linearity, this result needs only be proved for monomials / (z) = z" : it then reduces to 
identity (|II.5I) . □ 

We are now in position to prove the following results. 

Theorem 13. The equation UV.^i) has a solution 

hi = tho - [g, ho]„ (IV.7) 

n>0 

The corresponding value of the Hamiltonian h is 

^ = (IV.8) 

n>0 

Hence the coefficients An in jlV.^D o-nd Kn in UV.!^) are respectively equal to 

jo, n = 0 

n>0 

and 

An = 2 ^En ^ ^ 0 - 

Proof. Rewrite first the right-hand side of (IIV.7P as 

hi = tho — hoe~'^^ ■ 

deduce 

hi + = iho — -I- e'^ (i/iq — e“'^ 

= i{ho + e«/ioe-«) - z (e«^hoe-«^ + • 

Applying PV.6I) with z = 0 and the function f (S) = reduces the second term in the 

right-hand side to 2ho so that 

h-i -I- e‘^hie~‘^ = i (—ho + e'^hoe~‘^) 
and hi is a solution to PV.2I) . The corresponding value of h is 

h = ei {ho + ihi) e~i = (ho — ho + e'^^hoe”'^^) e~i 

Expanding the exponentials in the formula gives the result (IIV.8|) . □ 


V. CONCLUSION 

We have studied three identities on commutators and anti-commutators in which the Euler and 
Bernoulli numbers play a fundamental role. An open problem and desirable result at this point would 
be to link these identities with the combinatorial meaning of the Euler numbers: see 0 for a possible 
application. 
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